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Abstract 

We propose nonlinear integral equations for the finite volume one-particle 
energies in the 0(3) and 0(4) nonlinear fi-models. The equations are writ- 
ten in terms of a finite number of components and are therefore easier to 
solve numerically than the infinite component excited state TBA equations 
proposed earlier. Results of numerical calculations based on the nonlinear in- 
tegral equations and the excited state TBA equations agree within numerical 
precision. 



1 Introduction 



A better theoretical understanding of finite size (FS) effects is one of tlie most 
important problems in Quantum Field Theory (QFT). The study of FS effects is a 
useful method of analysing the structure of QFT models and it is an indispensable 
tool in the numerical simulation of lattice field theories. 

Finite size effects can be studied through the volume dependence of the mass 
gap of the theory, the usefulness of which is demonstrated P by the introduction 
of the Liischer-Weisz- Wolff (LWW) running coupling that enables the interpolation 
between the large volume (non-perturbative) and the small volume (perturbative) 
regions in both two-dimensional sigma models and QCD. 

The study of the Liischer-Weisz- Wolff running coupling is useful in the two- 
dimensional 0{N) nonlinear sigma (NLS) models, because according to recently 
performed high precision Monte Carlo measurements of the LWW running coupling 
[2j the cutoff effects look linear in these models, in contrast to perturbative consid- 
erations. The knowledge of the exact value of the LWW coupling enables one to 
make better fits for the cutoff effects, and thus to determine more accurately the 
functional form of the lattice artifacts. 

Our aim in this paper is to propose nonlinear integral equations (NLIEs) for 
the one-particle states of the 0(3) and 0(4) NLS models that allow for a fast and 
accurate numerical calculation of the LWW coupling. Although this calculation 
has recently been done in p[ using the excited state Thermodynamic Bethe Ansatz 
(TBA) technique jH El El 131310!, the difficulties corresponding to the infinite number 
of components in the TBA equations make the numerical calculations slow and at 
the same time restrict their accuracy too. Therefore it is desirable to work with the 
more convenient NLIE technique. 

Another point of the construction of these excited state NLIEs is to demonstrate 
that, similarly to the case of the sine-Gordon model JOl E], the NLIE technique 
can be extended to describe the finite size excited states also in the family of NLS 
models. 

The infinite set of TBA equations for the ground state of the 0(3) and 0(4) 
NLS models were given in ^21 ^| and ^1 1^ respectively. The derivation of the 
equations were based on the fact that NLS models can be represented as (limits of) 
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certain perturbed conformal field theories [121 CHI Cll CHI ■ 

The TBA description of the excited states is less systematic in the continuum 
models as it is for the ground state problem. Although a lot of different methods 
have been worked out to obtain excited state TBA equations in different models 
lU El El 13 IHl E! , a general construction has not discovered, yet. The generalization of 
the TBA equations of the 0(3) and 0(4) NLS models to one-particle excited states 
was proposed recently [Sj. The sigma- model TBA equations consist of infinitely 
many components, which makes their numerical analysis difficult. 

The NLIEs for the ground state of the 0(3) NLS model were proposed in |16j . 
based on the statement, according to which the 0(3) NLS model can be expessed as 
a certain limit of appropriately perturbed Zjv parafermion conformal field theories 

The ground state NLIEs for the 0(4) NLS model were derived in ^^1, using the 
integrable lattice regularization of the model [201 |2I1 1221 ■ 

Our main purpose in this paper to propose excited state NLIEs for the one- 
particle states of the 0(3) and 0(4) NLS models. This is achieved in sections 4 and 
5 using the assumption that the excited state NLIEs differ from the ground state 
ones only in additional source terms, and if it is necessary in additional quantization 
conditions. 

The paper is organised as follows. In section 2 we recall the TBA integral 
equations, Y-systems and NLIEs corresponding to the ground state problem. In 
section 3 we briefly summarize the one-particle TBA equations of the models. In 
sections 4 and 5 we propose excited state NLIEs for the 0(4) and 0(3) NLS models 
respectively. Numerical solutions of the NLIEs and their comparison to perturbation 
theory and to TBA results are discussed in section 6. Finally our conclusions are 
summarized in section 7. 

2 The ground state problem of the 0(4) and 0(3) 
NLS models (TBA and NLIE) 

In this section we give a short review of the ground state TBA equations, Y-systems, 
and nonlinear integral equations (NLIEs) for the 0(4) and 0(3) NLS models. 
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The ground state energy of a two-dimensional integrable model enclosed in a 
finite box with periodic boundary conditions can be determined by the solutions of 
the TBA integral equations prHJ. The TBA equations of the 0(4) and 0(3) NLS 
models can be encoded in infinite Dynkin- diagrams. (See fig. la and 16) 

Q 2 
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Figure 1: Dynkin-diagrams associated with the Y-systems of the a.) 0(4) 
and b.) 0(3) cr-models. 

The unknown functions ya are associated to nodes of the Dynkin-diagram and 
the TBA equations are of the form [T^OITlllTHl 



ya{x)=exp^SaiV{x) + Y,Iab{K*\ogYb){x)\, a =1,2,... (1) 

where 



6=1 



Ya{x) = 1 + ya{x), K{x) = ^^Qg^^TT^ ) '^i^) = cosh y-xj , (2) 

oo 

the * denotes the convolution, i.e. (/ * g){x) = J dy f{x — y)g{y), m is the mass 

— oo 

gap in infinite volume, / is the box size and lab is the incidence matrix of the Dynkin- 
diagram. The TBA equations of the 0(4) NLS model correspond to the diagram 
shown in figure la \l^. where the oriented double line at the beginning of the 
diagram means 

Ii2 = 2, hi = 1. (3) 

The TBA equations of the 0(3) NLS model are encoded into a D^o diagram shown 
in figure The ground state energy can be calculated from the solutions of the 
TBA equations O [13 El CSI 

oo 

^o(0 = -j I dx cosh (^) logri(x). (4) 

— oo 

The solutions of the TBA equations also satisfy the so called Y-system equations 

m 

ya{x + i) ya{x-i) =Y[Yb{xY^\ (5) 
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a{x) , a(x) , F (x) *i a{x), a{x), F{x) 

Figure 2: Graphical representation for the NLIEs of the a.) 0(4) and b.) 
0(3) cj-models. 

There is another alternative to calculate the ground state energy of the 0(4) NLS 
model, namely the nonlinear integral equation technique. The nonlinear integral 
equations for the ground state of the 0(4) NLS were derived in [T2j using the light- 
cone lattice approach of the model [20]. The nonlinear integral equations in this 
case are of the form 

\ogyi{x) = V{x) + 2{K+^ *\ogU){x) + 2{K-'' *hgU){x), 

log a{x) = {F* log U) (x) - (F+2(^-^) * log U) (x) 

+ {K-U\ogYi){x), (6) 

loga(x) = (F*logf/)(x) - (F-2(i-T) *logf/)(x) 

+ {K+^U\ogY^){x), 

U{x) = l + a{x), U (x) = 1 + a{x) , Yi{x) = 1 + yi{x) , 



dq e-t'^l-^"- 



where 0<7<l/2isan arbitrary fixed parameter, 

dq e-'9l-^' 
2tt 2 cosh(g) 

— oo 

and we have used the notation 

f^\x) = fix±tv). (8) 

In this case the form of the ground state energy is the same as it is in the TBA case 
(HI). Equations © contain only three real unknown functions because yi{x) and 
Yi{x) are real and a{x) is the complex conjugate of a{x), therefore they serve an 
efficient basis for numerical calculations. The graphical notation of these equations 
is depicted in figure Eh- The big "bubble" denotes the complex auxiliary functions 
which resum the contributions of those TBA nodes, which are inside it. In our 
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notation the names of the complex unknown functions and the kernel function are 
indicated. 

In the eqs. (jH)) were derived on a Bethe Ansatz solvable lattice and certainly 
the TBA equations of the model can be also derived from the Bethe Ansatz solution 
of the model. Thus, it turns out that the function yi{x) of eqs. (jH)) are exactly the 
same as the function yi{x) of the TBA equations (Q). Furthermore the connection 
between the complex a{x) and a(x) variables and the TBA variables can be expressed 
by the following formula 

U{x + i-f) U{x-i-f) =Y2{x), (9) 

which is a very important formula from the NLIE technique point of view, because 
this relation allows one to reduce the infinite component TBA system to a finite 
component NLIEs. 

NLIEs are also available for the ground state of the 0(3) NLS model. In ^H] 
NLIEs were proposed to a class of perturbed parafermion conformal field theories, 
which reduce to the 0(3) NLS model in a certain limit yL2j. Taking the appropriate 
limit the NLIEs take the form 

log yi (x) = V{x) + log y2 (x) , 
logy2(x) = *\ogU){x) + {K-'' *\ogU){x), 

log a{x) = {F* log U) (x) - (F+2(i"7) * log U) (x) 

+ {K~^ *\ogYi){x) + {K-^ *\ogY2){x), (10) 

loga(x) = (F*log[/)(x) - (F'2(^-^Ulogf/)(x) 

+ {K+^ *\ogYi){x) + {K^^ *\ogY2){x), 

U{x) = l + a(x), f7(x) = 1 + a(x), Ya{x) = 1 + ya{x), a = 1,2 

where 0<7<l/2isan arbitrary fixed parameter, the kernel and source functions 
are the same as in (jS)),© and the ground state energy is given by the formula (jl} as 
in the TBA case. The graphical representation of these equations is given in figure 

Although these equations are only conjectured ones, one can recognise that the 
functions yi{x) and y2{x) are the same as the ones of the corresponding TBA equa- 
tions (^, a{x) and a{x) are the complex conjugates of each other and the connection 
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between the TBA and NLIE variables can be expressed (similarly to the 0(4) case) 
by the formula 

U{x + i-f) U{x-i-f) = Ys{x), (11) 
which somehow ensures the reduction of the infinite component Y-system. 

3 TBA equations for the one-particle states of the 
0(4) and 0(3) NLS models 

In this section we briefiy review the first excited states TBA equations of the 0(4) 
and 0(3) NLS models proposed in [Sj. 

In jSj the main assumption was (based on previous experience with the sine- 
Gordon case) that the Y-system (jSj) describes not only the ground state of the 
model, but remains valid also for the first excited state. Then Liischer's mass gap 
formula valid for asymptotically large volumes, was used to determine the 
infinite volume solution of the excited state Y-system, which is sufficient to derive the 
one-particle excited state TBA equations. It is important from the NLIE technique 
point of view that the infinite volume solutions of the Y-system can be expressed 
by an infinite volume t-system, which is of the form 

tp{x + i) tp{x - i) = B{x - ip) B{x + ip) + tp^i{x) t.p^i{x), p = l,2, ... (12) 

where 

to{x) = 0, (13) 

B{x) = X, in the 0(4) case, and (14) 

B{x) = x^ + l, in the 0(3) case. (15) 

The solution of (fT^ can be represented in the following form 

p 

tp(a;) = ^Af^(x), Xl^^x) = B [x + i{2l - p - 1)] . (16) 
1=1 

Using the infinite volume t-system ()12|1 . the infinite volume solutions of the Y-system 
can be given as 

^oo(^) ^ tp^i{x)tp+i{x) Y^(x)= ^p(^ - ')^p(^ + p=l,2,... 

^ B{x — ip)B{x + ipY ^ B{x — ip)B(x + ip)^ ' ' 

(17) 
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The infinite volume Y-system of tfie 0(3) and 0(4) NLS models [H] are of tlie form 

y^{x + i) y^{x -i) = Y^_,{x) Y^^,{x) p = 2, 3, ... (18) 

where y^{x) = 0, because the Y-system element corresponding to the massive node 
of the Dynkin-diagram tends to zero in the infinite volume limit. The usefulness of 
these infinite volume solutions is that, they give the qualitative position of the zeroes 
of the ya{x) functions of the Y-system, which is sufficient jHl for the derivation 
of excited state the TBA integral equations. Thus the one particle TBA equations 
of the 0(4) NLS model are of the form 

yaix) = T^ix)expi^SaiVix) + J2lab{K*logYh)ix)\ , a =1,2,... (19) 

where 



6=1 



t{x) = tanh (^^j , (20) 
and lab the incidence matrix of figure la. The energy of this state takes the form 

oo 

in f / 7TX \ 

Ei(/) = m - — / dxcosh (^— j logyi(a;). (21) 

— oo 

Similarly the one particle TBA equations of the 0(3) NLS model can be read off 
from ^ and take the form 

y,ix) = e^(^)y2(x), 

y^ix) = T{x-h) r{x + h) exp{{K*\ogY3){x)} (22) 

I E ^ab iK*logYh)ix)\ 

ya{x) = t{x - ha-l) t{x + ha-l) t{x - /l^+l) t{x + /l^+i) 6 '■f=i J 

a = 3,4,... 

and in addition to these integral equations one has to impose quantization conditions 
for the ha zeroes of the ya{x) functions to ensure that ya{^ha ± i) = —1, which is 
the consequence of the (jSI) Y-system relations. These quantization equations can be 
obtained by the analytic continuation of (j22|) and are of the form 



^ oo 

'j{hs-hs-i)+'jihs+hs-i)+'j{hs-hs+i)+'j{hs+hs+i) — : V'/st iK~^*\ogYb){x) = vr, 

(23) 
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for s = 3,4, ... where h2 = 0, 7(x) = 2arctanr(x) and 



oo 



{K-U\ogYb){x) = iV 



I 



dy logYb{y) 
4 sinh |(a; — ' 



(24) 



— oo 



is a principal value integration. 

Our main goal in this paper to propose NLIEs equivalent to these infinite com- 
ponent excited state TBA equations. 



4 Nonlinear integral equations for the first ex- 
cited state of the 0(4) NLS model 



In this section we propose nonlinear integral equations for the first excited state 
of the 0(4) NLS model. The proposal is based on the assumption that the excited 
states NLIEs differ from the ground state ones only in additional source terms, plus if 
it is necessary in quantization conditions. In addition we assume that the functional 
form of these additional source terms does not depend directly on the volume. They 
can depend on the volume only through some complex objects, on which we impose 
quantization conditions. Such objects may be the zeroes of the Y-system elements, 
or some roots of the Bethe Ansatz equations as it happens in the sine-Gordon case 
jllj . These assumptions are valid for the NLIEs in the sine-Gordon model, where 
the form of the excited state NLIEs of ref. fl] differ from the ground state NLIE 
of ref. ^0] only in additional source terms and in quantization conditions. 

In this model one can even give a physical interpretation of the various terms of 
the equation. The integral term represents the contribution of the infinite number 
of real Bethe Ansatz roots. The presence of these roots can be regarded as filling the 
Dirac sea and thus the presence of this term in the NLIE of the sine-Gordon model 
is somehow related to the vacuum structure of the theory. For an excited state of 
the model one has to consider a slightly modified configuration of the Bethe Ansatz 
roots. It turns out that in the sine-Gordon model these modifications can be taken 
into account by adding source terms to the NLIE of the ground state problem [TT] . 
Similarly, the excited state TBA equations of the sine-Gordon model differ from the 
ground state ones only by some source terms ^Ij. Finally, additional support for 
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our assumptions comes from a direct calculation in the higher-spin vertex model, 
which reduce to the 0(4) NLS model in the infinite spin limit. These calculations 
show that excited state NLIEs differ from the ground state ones again only by some 
source terms, which depend on the volume only through some complex objects, on 
which quantization conditions have to be imposed. (As an additional example see 
PH] for the S=l special case.) However, we cannot directly derive the one-particle 
NLIEs of the 0(4) NLS model from the (limit of the) higher-spin vertex model, since 
the solvable lattice model is describing only a subsector of the Hilbert space of the 
model [201 1^ and unfortunately the one-particle states are not in this subsector. 

Accepting these assumptions our task is to find out the necessary source terms 
of the excited state NLIEs. This is achieved similarly to the case of the derivation of 
the excited state TEA equations of the model, where the infinite volume solutions of 
the equations allowed one to deduce the TEA integral equations . The unknown 
functions of the excited state problem will be denoted in the same way as in eqs. 
©. Assuming that the relation Q holds also for the first excited states, we can use 
this relation in the equation for the massive mode in (^. This leads to the following 
Ansatz for the first excited state NLIE problem: 

logyi(x) = V{x) + \ogT^{x) + 2{K+^ *\ogU){x) +2{K-^ *\ogU){x), 

loga(x) = J^a{x) + {F*\ogU){x)~{F+'^^^-^^ *hgU){x) 

+ {K-^*logY^){x), (25) 

loga(x) = J^^j(x) + (F*logf/)(x) - (F-2(i-T) *logf/)(a;) 

+ {K+^*\ogY,){x), 

U{x) = l + a{x), U{x) = l + a{x), Yi{x) = 1 + yi{x) , 

where 0<7< 1/2 is an arbitrary fixed parameter, the energy formula identical to 
()21|) and our task is to determine the presently unknown J^a{x) and J^a{x) source 
functions. Using the assumption that these functions do not depend on the volume 
directly, one can determine them from the infinite volume solution of the proposed 
equations, which can be determined from the solution of the infinite volume t-system 
()12|1 and Y-system (fTH|) . Using the construction of pT| 1^ . the infinite volume 
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solutions of the Ansatz take the form 

^ Xf\x + i-i^) ^ ti{x-i-f) ^ x-i'y 
''^^''^ X?\x + ^-^^) 5[x + .(2-7)] x + ^(2-7)' ^ 

_ Xfjx-i + i-f) ^ ti(x + i7) ^ a: + z7 

''"^^''^ Af)(a:-. + z7) i?[^-^(2-7)] x-.(2-7)' ^ 

U^{x) = 1 + a^{x) = 2 "^^'i^'^l , (28) 

X + t{2 - 7) 

f/oo(x) = l + aoo(x) =2 (29) 

X — 2(2 — 7) 

and y^{x) = 0. One can check that a{x) is the complex conjugate of a{x) and 
that these solutions satisfy the relation Uoo{x + 27) Uoo{x — 27) = Y^{x), which 
is the infinite volume limit of (jHl). Once one knows the infinite volume solutions 
of the Ansatz one can substitute these formulae into them, and can compute 
explicitly the unknown source functions. The most convenient way of this is to take 
the derivative of eqs. (f^K|) and work in Fourier space, and in the end return to the 
coordinate space and integrate once. After these simple calculations one gets 

J^a{x) = x{x + t{l - 7)) + Xk{x + t{l - 7)), (30) 



-^a(x) = -x{x - i{l - 7)) - Xk{x - i{l - 7)), (31) 



where 



X 00 

-kl 



X{x) = 2ni I dy F{y) = i [ — sm{qx) ^ . (32) 
J J q 2cosh(g) 

-00 

X 

Xk{x) = 2,711 j dy K{y) = i arctansinh (^-^^ ■ (33) 







Note that these functions are nothing but the odd primitive of the kernel functions of 
the integral terms of the ground state equations (jHI), which is similar to what happens 
for the NLIE of the sine-Gordon model in the case of the one hole excitation [llj. 

We solved these NLIEs numerically for a number of cases and the results we 
found always agreed with those obtained from the numerical solution of the TBA 
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equations p9|) (see section 6.). This means that the two methods are equivalent. 
On the other hand the TBA equations were tested previously by using Liischer's 
formula, Monte Carlo measurements P| and 3- loop perturbation theory Thus 
we are confident that eqs. ()25|) with source terms 1)301311) correctly describe the 
one-particle excited states of the 0(4) NLS model. They are superior to the TBA 
equations since they contain only three real unknown functions leading to faster 
convergence in numerical calculations. 

5 Excited states nonlinear integral equations for 
the 0(3) NLS model 

In this section following the method described in the previous section we propose 
NLIEs for the first excited state of the 0(3) NLS model. According to our assump- 
tions the Ansatz for the NLIEs is of the form 

\ogyi{x) = V{x) +\ogy2ix), 

log 7/2(0;) = logT{x - hs) + logT{x + hs) + {K'^'^ *\ogU){x) + {K~'^ *logU){x), 

\oga{x) = J'a(x) + (F*logf/)(x) - (F+2(i-7) *logf7)(x) 

+ {K'^ *\ogYi){x) + {K-^ *\ogY2){x), (34) 

log a{x) = J^-a{x) + {F* log U) (x) - {F-^C^-^'^ * log U) (x) 

+ (K+^*logFi)(x) + (ir+^*logF2)(x), 

U{x) = l + a{x), U{x) = l + a{x), Ya{x) = 1 + ya{x), a = 1,2 

where 0<7<l/2isan arbitrary fixed parameter, yi{x) and y2{x) are the same 
variables as the ones in the excited state TBA eqs. ((221), and is the zero of 
the ysi^x) Y-system element of the excited state problem (i.e. y^i^ih^ ± i) = — 1), 
therefore an additional quantization equation must be imposed on this zero. This 
quantization condition can also be derived from the infinite volume solution of our 
equations, and will be discussed at the end of this section. The energy expression 
in turn is the same as the TBA one (PT|) . 

According to our method first we have to determine the infinite volume solution 
of our Ansatz This can be achieved using the solutions of the infinite vol- 



11 



ume t-system p2lll6|) . and applying the construction of refs. PT| I25j. After some 
straightforward calculations one gets 

xf\x + i — 27) + X^2^ (x + 2 — 27) (x — 27)^ 



Xf\x + ix + ^(4 - 7)) ■ (x + 2(2 - 7)) ' 

(35) 



/ \ -^2^^ (x - i + + A?^ (x - z + 27) ^ (x + 27) 



Xf\x + (^ - ^(4 - 7)) ■ {x - i{2 - 7)) ' 

(36) 

U {x)-l + a ix) 3 (^ + ^o + ^(l-7))-(x-/^o + ^(l-7)) . 
f/.(x) ^ 1 + a.(x) ^ 3 + 

(x -z(4-7)) ■ (x -2(2 -7)) 

00/ N ^ ^i(a^) ^3(3^) ^ o (a^ + /^o) {x - /iq) 
^2 5(x + 22) 5(x - 22) (x + 32) (x - 32) ' ^ ' 

Y^ix) - 1 + V^ix) - ^2(x + 2) t2(x-2) _ (a: + 2) (x - 2) 

^2 (Xj - i + 2/2 (X) - ^ _ 2z) - ^ (x + 32) (X - 32) ' ^^^^ 



where = lim /i3 = ^^5/3, the infinite volume limit of the zero /i3, and certainly 
; / / ■ = in the infinite volume limit. The infinite volume solutions ()35ll4()j) also 
satisfy the important relation f/oo(x + i'y)Uoo{x — 27) = Y^{x), which allows one to 
cut the infinite TBA equations. 

Having the infinite volume solutions of our Ansatz, using the method detailed 
in the previous section one can calculate the unknown source terms (at least their 
infinite volume limit). After simple calculations one gets 

:Fa{x) = x{x + ho + i{l- 7)) + x{x-ho + i{l- 7)) + 2xk {x + 2(1 - 7)) , (41) 



Mx) = -x{x + ho- 2(1 -j))-x{x-ho- 2(1 - 7)) - 2xK {x - 2(1 - 7)) , (42) 

according to our assumption these source terms may depend on the volume through 
the objects, which can be found in their arguments, namely through the volume 
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dependence of Hq. Therefore we replace by its finite volume value in ()41I42|) . 
and according to our conjecture that will be the form of the source functions for 
finite volume, 

^a{x) = x{x + h + i{l- 7)) + x{x-h + i{l- 7)) + 2xk {x + 1(1 - 7)) , (43) 



^a(x) = -x{x + h- 2(1 -l))-x{x-h- 2(1 - 7)) - 2Xk {x - 2(1 - 7)) • (44) 

The eqs. (jH^ must be supplemented by the quantization condition of the zero h^, 
similarly as it was in the TBA case. The quantization condition can be found out 
from the infinite volume solutions ()35II4()|1 . One can see from these solutions that 
f/oo(±/?.o ^ ^(1 ^ 7)) = 0. Assuming that this relation remains valid also for finite 
volume, one can infer that f/(±/i3 — i{l — 7)) = 0, from which follows that 

a(/i3-z(l-7)) = -l- (45) 

Taking the logarithm of this, and using eqs. with ()43I44|) . one gets the quanti- 
zation condition for 

Aih) = vr, (46) 

where 

1 1 2 1 

A{x) = - xix-h) + - x{x + h3) + - xk{x) + - iK-UlogYi){x) (47) 

2 I t I 

+ - (K-^ * logF2)(a;) + - (F^^^^^) * \ogU)(x) - - (F+^i-^) * logm(x). 
22 2 

To summarize the eqs. (jH^ with the source functions ()43l44j) . and supplemented by 
the quantization conditions ()4fip47j) make up the one-particle NLIEs of the 0(3) NLS 
model. In these equations one can also see that the source functions are nothing but 
the odd primitives of the kernel functions, as in the 0(4) case. Nevertheless in this 
case the situation is a bit more complicated, because they contain some zeroes in 
their argument, namely the function x{x) contains the zero /13, which corresponds to 
the zero of the 2/3(0;) (in the y3(±/i3±z) = — 1 sense), and the function Xk{x) contains 
the zero /12 = 0, which corresponds to the zero of the 2/2(2;) (in the 2/2(^^2 ±2) = — 1 
sense), but for this object we do not need to impose quantization equation because 
of symmetry reasons. 
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Following from the construction of these equations, one can analytically verify by 
solving iteratively these equations for large / that they give back correctly the same 
leading order correction to the infinite volume mass gap as it is predicted by the 
Liischer's formula [22] ■ For not very large volumes we tested these NLIEs through 
numerical calculations, and we found that these equations give the same numerical 
results as the earlier proposed excited state TBA equations fj22ll24|) (see section 6.). 
Because the corresponding TBA system was tested by the results of 3-loop order 
perturbation theory and lattice Monte Carlo measurements 0, thus our new 
excited state NLIEs are also consistent with these methods, and correctly describe 
the first excited state energy of the 0(3) NLS model. 

6 Numerical results 

In this section we perform numerical checks on our conjectured one-particle NLIEs. 
We solved numerically both the TBA equations and the NLIEs for the ground state 
and for the one-particle states at some values of the volume, and we compared the 
numerical results of the two different methods. Moreover we solved our conjectured 
NLIEs in the deep ultraviolet region so as to be able to compare our numerical 
results to the predictions of the asymptotically free perturbation theory. 

The numerical method used for solving the TBA equations is described in 
The numerical solution of the corresponding NLIEs is rather similar. Namely, we 
solve numerically the NLIEs by iteration. The starting point for the iteration is the 
I — > oo solution of the equations and the procedure converges rapidly. The one- 
particle state problem for the 0(3) NLS model is more involved since here one step 
of the iteration includes the calculation of the integrals occurring in ()34|) together 
with the calculation of the zero (h^) from the quantization condition ()4(ip47|l . Again, 
the starting point of the iteration procedure is given by the I — > oo solution, both 
for the unknown functions and for the position of the zero h^. 

We used Simpson's formula, and a cutoff A for the numerical evaluations of the 
integrals running from —oo to +oo, in such a way, that in the region of |x| > A, the 
unknown functions are approximated by their infinite volume limit. The magnitude 
of the cutoff A at the values of the volume under consideration has to be 100 in the 
0(4) case and 150 in the 0(3) case, so that the 9-digit numerical precision could be 
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ml 


Eo (NLIE) 


^0 (TBA) 


2 


-0.1620289681(1) 


-0.16202897(1) 


1 


-0.6437745719(1) 


-0.6437746(1) 


10-1 


-11.273364587(1) 


-11.273364(1) 


10-2 


-127.22634373(1) 


-127.2263(1) 


10-3 


-1343.4090793(1) 


-1343.409(1) 


10-^ 


-13865.238816(1) 




10-^ 


-141563.8217(1) 




10-^ 


-1436683.423(1) 





Table 1: NLIE and TBA results for the ground state energy in the 0(4) NLS model. 



ml 


Eo (NLIE) 


Eo (TBA) 


2 


-0.1228466915(1) 


-0.1228466(1) 


1 


-0.4862495672(1) 


-0.4862496(1) 


10-1 


-8.006985662(1) 


-8.006985(1) 


10-2 


-87.63570019(1) 


-87.6357(1) 


10-3 


-913.9547387(1) 


-913.954(1) 


10-4 


-9374.188294(1) 





Table 2: NLIE and TBA results for the ground state energy in the 0(3) NLS model. 

reached. Our numerical results are summarized in Tables 1-4. 

One can see from these numerical data, that the numerical results served by our 
NLIEs agree with those of the TBA equations within the numerical precision. 

Now we are able to compare our numerical results to those of the asymptotically 
free perturbation theory. Having the numerical values of the ground state energies 
and the one-particle state energies, we can compute numerically the dimensionless 
finite volume mass gap (LWW coupling) 

z{ml) = I ■ [E,{1) - Eo{l)], (48) 

for which perturbative results are also available up to 3-loop order [26^ . The pertur- 
bative formulae neccesary to compute the 3-loop perturbative mass gap can be found 
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ml 


El (NLIE) 


El (TBA) 


2 


0.9923340593(1) 


0.99233406(1) 


1 


0.9383970591(1) 


0.93839706(1) 


10-^ 


-3.004108884(1) 


-3.0041089(1) 


10-2 


-69.83802786(1) 


-69.838028(1) 


10-^ 


-901.2815867(1) 


-901.28159(1) 


10-^ 


-10260.214298(1) 




10-5 


-111091.0324(1) 




lO"'^ 


-1172575.496(1) 





Table 3: NLIE and TBA results for the one-particle state energy in the 0(4) NLS model. 



ml 


El (NLIE) 


El (TBA) 


2 


1.02169721(1) 


1.0216972(1) 


1 


1.084208673(1) 


1.084208(1) 


10-1 


0.77721084(1) 


0.77718(1) 


10-2 


-23.6407101(1) 


-23.643(1) 


10-^ 


-406.195912(1) 


-406.23(1) 


10-^ 


-5150.21619(1) 





Table 4: NLIE and TBA results for the one-particle state energy in the 0(3) NLS model. 
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ml 


z{l) (NLIE) 


3-loop PT 


2-loop PT 


10-1 


0.826925570(1) 


0.826130 


0.8252260 


10-2 


0.573883159(1) 


0.5737662 


0.5735488 


10-3 


0.442127493(1) 


0.4420969 


0.4420193 


10-4 


0.360502452(1) 


0.3604916 


0.3604571 


10-5 


0.30472789(1) 


0.3047233 


0.3047056 


10-6 


0.26410793(1) 


0.2641057 


0.2640957 



Table 5: NLIE and PT results for the finite volume mass gap z(mZ) in the 0(4) NLS model. 



ml 


z{l) (NLIE) 


3-loop PT 


2-loop PT 


10-1 


0.878419650(1) 


0.876058 


0.873458 


10-2 


0.639949901(1) 


0.639645 


0.638874 


10-3 


0.507758827(1) 


0.507669 


0.507358 


10-4 


0.422397210(1) 


0.422363 


0.422212 



Table 6: NLIE and PT results for the finite volume mass gap z{ml) in the 0(3) NLS model. 

in The comparison of our numerical results to the predictions of the perturba- 
tion theory can be found in Tables 5 and 6. We listed both the 2-loop and 3-loop 
perturbative results so that one can infer to the accuracy of the 3-loop perturbation 
theory at the scales under investigation. From this comparison, at very small ml 
values one can experience a very nice 4-digit agreement in the 0(3) case, and an 
almost 6-digit agreement in the 0(4) case, which is non-trivial, since the conjecture 
of our NLIEs was based on the large volume asymptotics of the unknown functions. 
This perfect agreement makes us confident that our conjectured one-particle NLIEs 
describe the exact finite volume one-particle energies in the 0(3) and 0(4) NLS 
models. 
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7 Summary and Conclusions 



In this paper we proposed NLIEs for the one-particle states in the 0(3) and 0(4) 
NLS models. The form of these excited state NLIEs are based on the assumption 
that they differ from the ground state ones only by some source terms, which may 
depend on the volume through some objects in their argument, on which extra 
quantization conditions must be imposed. This assumption is mainly motivated 
by the form of the NLIEs in the sine-Gordon model and by direct calculations in 
the higher-spin vertex model (which go to the 0(4) NLS model in the infinite spin 
limit). Accepting these assumptions and starting from the explicit infinite volume 
solution of the first excited state Y-system of the models we were able to find the 
infinite volume limit of the conjectured equations. This is sufficient to determine 
the source terms and the quantization conditions. It is interesting to note that in all 
cases the source terms are the odd primitives of the kernels occurring in the integral 
terms of the equations, just as in the case of the NLIE in the sine-Gordon model. 
We have solved these equations numerically and found that results agree with those 
obtained previously from numerical solution of the excited state TBA equations and 
in the deep ultraviolet regime the numeical results also agree with the predictions 
of the 3-loop perturbation theory. This agreement is convincing evidence for the 
correctness of our integral equations. The advantage of using these NLIEs (instead 
of the infinite set of TBA equations) is that here the number of unknown function 
is finite (and small). 

An interesting generalization of our results would be to propose excited state 
NLIEs for the (p(id,id,adj) perturbation of the SU{2)l x SU{2)k / SU{2)l^k models, 
which were investigated at rational level K in [TE] . 

It would also be interesting to extend the NLIE technique for all the excited states 
of the 0(3) and 0(4) NLS models and for such more complicated TBA systems which 
can be encoded into the products of two Dynkin- diagrams P7] . 
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